Nisho-1-2012 



Quark Pair Production in Expanding Glasma 

Aiichi Iwazaki 

International Economics and Politics, Nishogakusha University, 
6-16 3-bantyo Tiyoda Tokyo 102-8336, Japan. 
(Dated: September 4, 2011. Revised: January 13 2012) 

Glasma in high energy heavy ion coUisions is longitudinal classical color electric and magnetic 
' fields. The color electric field has been shown to produces quark and anti-quark pairs by the 

I , Schwinger mechanism and to oscillate with time in non-expanding glasma, that is, plasma oscillation. 

■ On the other hand, in the expanding glasma we show that the field decreases with the plasma 

' oscillation. We can explicitly obtain the solutions representing such temporal behaviors in the 

system with r and rj coordinates. We show these results by using massless QED as a simplified 
model of QCD. 
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r-| . I. INTRODUCTION 

^ • Initial states of color gauge fields (glasma) produced immediately after high energy heavy-ion collisions have recently 
1^ ' received much attention. The gauge fields are longitudinal classical color electric E and magnetic B fields. The 
■ presence of such classical gauge fields has been discussed on the basis of a fairly reliable effective theory of QCD at 
high energy, that is, a model of color glass condensate (CGC)[ll, The glasma is generated by small x gluons in 
^■f^ , nuclei. It is expected that the decay of the glasma leads to thermalized quark gluon plasma (QGP). The study of 
^ ■ the decay of the glasma is very significant because the classical and quantum behaviors of the coherent non-Abelian 
On gauge fields have not been experimentally known. 

CO The glasma is homogeneous in the longitudinal direction but inhomogeneous in the transverse directions. Hence, we 

' may view that it forms electric and magnetic flux tubes extending in the longitudinal direction. Their field strengths 
[ I are given by the saturation momentum Qs in high energy heavy ion collisions such that eB ^ and eE ~ Q^; e(> 0) 
On denotes the coupling constant. In the previous papers[3|-|5| we have shown a possibility that the famous Nielsen- Olesen 
, instability |[6| makes the color magnetic field B decay. The possibility has been partially confirmed by the comparison 
between our results and numerical simulations 0, Such decay is a first step toward the generation of QGP. ( The 
gauge fields produced initially in the model of color glass condensate have very random spatial distribution without 
no typical length scales[9]. The Nielsen- Olesen unstable modes do not arise under such random configulations of the 
fields. But soon after the production they would rapidly evolve to form a smooth distribution with the typical length 
scale Qj^ owing to the nonlinear interactions among the fields O. Then, the unstable modes can arise under the 
smooth background gauge fields. ) 

Furthermore, we have also discussed pll - [l3| the decay of the color electric field; the decay is caused by Schwinger 



mechanism IJ] , that is, the pair production of quarks and anti-quarks. The mechanism has been extensively 
explored [l5 since the discovery of the Klein paradox. A new feature in the glasma is that it is composed of not 
only electric field but also magnetic field collinear to the electric field. Under such a circumstance there are only 
several studies of the mechanism involving the back reaction of the produced particles on the electric fieldjl2l. [l5lll6j|. 
It has been found that the number density of quarks and the electric field oscillate in time, namely the plasma os- 
cillation when the quarks are free after their production. On the other hand, it has recently shown TtII that when 
the quarks interact with each other in heat bath, the electric field decays without the oscillation, while their number 
density increases and is saturated. These results have been obtained in a system of non-expanding glasma using QED 
for simplicity. 

In this paper we analyze the behaviors of the electric field and the quark number density in the expanding glasma. 
We only consider the case of non-interacting quarks after their production. The pair production in the expanding 
glasma has already been discussed fl 8], but there is still no analysis in the case that both electric and magnetic 
fields are present. In order to describe the expanding glasma, we choose a coordinate system such as r = yjt"^ — 
and r\ = log ^Jit -\- z)/{t — z) instead of t and z in Cartesian coordinate. ( r is a proper time associated with a 
particle whose longitudinal coordinate is specified by ry = constant. ) For simplicity, we analyze U(l) gauge theory 
instead of SU(3) gauge theory. Thus, we use the terminology of electrons and positrons. They are assumed to be 
massless particles. This is because the corresponding masses of quarks can be neglected in the glasma. The typical 
scales in the glasma is given by the saturation momentum Qs, which is much larger than the mass of the quarks ; 
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Qs ~ IGeV ~ 2GeV for RICH or LHC. We also assume that only relevant states in the pair production are those 
in the lowest Landau level with their energies given by \pz\ in the Cartesian coordinate where Pz is the momentum 
parallel to B. The assumption holds only in the limit of eB oo, although both eB and eE are of the order of in 
the real glasma. We discuss that the assumption is approximately correct even when eB and eE are the same order 
of magnitudes as each other. We find that both the electron number density and the electric field decrease oscillating 
with time. The decrease is caused by the expansion of glasma. 

In order to show these results, we use the chiral anomaly. The chiral anomaly is a quite powerful tool for the 
discussion of the Schwinger mechanism. It is effective only when both coUinear electric and magnetic fields are 
present. When the magnetic field is very strong, the produced particles occupy only the states in the lowest Landau 
level. It implies that the transverse motion of the particles are frozen and only the longitudinal motion is allowed. 
Thus, the particle production effectively takes place in two dimensional space-time. Then, the form of the chiral 
anomaly is much simplified. The pair production can be discussed without detail calculations but simply by solving 
the anomaly equation. 

In the next section [11] we explain how the pair production effectively arises in two dimensional space-time when the 
magnetic field is very strong. We define a two dimensional field operator using creation and annihilation operators of 
the states in the lowest Landau level. Electric charge and chiral current of the states are represented with the field 
operator. The analysis in this section is performed in the Cartesian coordinate. In the section [TTTl we show the utility 
of the chiral anomaly for the pair production by using the field operator. In the section IIVI we apply our method 
to analyse the pair production in r and 77 coordinates. We find that the electric field decreases oscillating with time 
owing to the back reaction and the expansion of the glasma. We summarize our results in the section [V] 



II. THE LOWEST LANDAU LEVEL AND REDUCTION OF SPATIAL DIMENSIONS 

First, we explain solutions of Dirac equation under the homogeneous magnetic field in the Cartesian coordinate and 
the reduction of the four dimensional operators to the two dimensional ones. We use the gauge potential A — (0, xB, 0) 
of the magnetic field B — {0,0, B). Then, the solutions representing massless fermions in the lowest Landau level 
which we are only concerned with, are given in the following. 



"if = Noexp{~iEpt + ipyy + ipzz)exp(- — {x ^)^)u{p) with u{p) 



where Ep = i\pz\ represents the energy and iVo denotes a normalization constant. On the other hand, the states in 
higher Landau levels have the energies ±-\/ enB + pi where the integer n > 1 denotes Landau level. Thus, in the limit 
ol eB — >■ 00, only relevant states in the pair production are the states with the energies \pz \ in the lowest Landau level. 
( In addition to the states there are other states in the lowest Landau level whose energies are given by ^ eB + pi . 
The states are those with magnetic moments anti-parallel to B, carrying much higher energies than the energies of 
the above states. These states are also irrelevant in the limit. ) Hereafter we only consider the states in the lowest 
Landau level, which are the most relevant states in the Schwinger mechanism when B ^ E. 

The transverse motions are described by the factor exp^+ipyu) exp(— ei?(x — Py/eB)'^/2) of the wave functions in 
the lowest Landau level. The states are degenerate in the momentum py. The states in higher Landau levels are 
described by different functions of the transverse coordinates x and y. In the sense, the transverse motions are frozen, 
as long as we are concerned with such states in the lowest Landau level. Only dynamical behaviors of the physical 
quantities are allowed in the longitudinal direction. Namely, the states in the lowest Landau level are specified with 
quantum numbers pz and py , whose energies are given by \pz \ . Since the electric field parallel to B accelerates the 
particles, only the momentum pz increases or decreases, but the momentum py does not change. Thus, the mometum 
Py is conserved and the particles do not jump into higher Landau levels as long as their energies \pz\ < eB. ( Although 
the momentum pz increases by the acceleration of the electric field, the electric field decays with the back reaction 
until Pz becomes sufficiently large such as \pz\ > eB. ) Namely, the transverse motions specified by the coordinates 
X and y are frozen. Therefore, we only need to specify the quantum number pz in order to describe the states in the 
lowest Landau level. 

Since they are trivial in the transverse directions, we take an average over the transverse space of physical quantities, 
e.g. with gamma matrices O. Then, we find that the average over the transverse space of are given by 

ipOilj where the field operator ip is defined as. 




(1) 
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= y ^ ( exp{—iujpt + ipz)u{p)Ap + exp(iajp — ipz)u{p)B^^ (2) 

with Up = \p\ and p = pz, where Ap and Bp denote annihilation operators of particles and anti-particles with the 
momentum p, respectively, and satisfy the following commutation relations, 

ApA^j + A^jAp = {Ap, AJ} — {Bp, B^} — 2n5{p — q), the other commutators vanish. (3) 
Then, it follows that 

{i,\t,z),i;{t,z')} = 5{z-z'). (4) 
The field i}} satisfies effectively 2 dimensional Dirac equation, {ij'^dt + 'y^dz)ip = where 

1 

1 
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7" = n _1 > 7^ = , and 7^ = 



Note that the spinor and the gamma matrices are not two dimensional one. The term of "two dimensional" means 
that the field depends only on the coordinates t and z. When homogeneous electric field E — —dtA3{t) is present 
in addition to the magnetic field, the operators Ap and Bp depend on t, but keep the above commutation relations. 
They satisfy the equations, idtAp = eAsApp/lp] and idtBjj = eAzB'^p/\p\. 

Using the field il), we find that the expectation values of the electric charge e'ip^'^ip, electric current chiral 
charge ipj'^^^ip and chiral current ipj'^^^ijj are given respectively by 



/dp — P dip p 
— {np-hp), and = {e^pj^ip) = e {rip - hp) (5) 
ivr J Ztt LOp 

J5 ^ (V;7VV'>= / ^ — {np-hp) and J^^^ ^ (V^tV^) = / ^{up-hp), 
J Zn Up J Ztt 

where operator products are defined as normal ordered product with respect to Ap and Bp. We have taken the 
expectation values of these quantities by assuming appropriate states | ) relevant to the particle production under the 
electric field, which satisfy 



(AlAg) = 2ttS{p - q)np{t), (bIb,) = 2ttS{p - q)hp{t) and (^tfit) ^ (^SpAg) = (6) 

Obviously Ap ( Bp ) represents annihilation operators of positrons ( electrons ) with the electric charge e > ( 
— e < ); np{t) represents momentum distribution of the positrons, while hp{t) represents that of electrons. Both of 
them depend on time t, in general. 

We denote the number density of positrons by iVp = J |f and of electrons by A^e = / ^tt-p respectively. Then, 
we find that 

jO = eNp - eNe, J^'^ ^ Np - N^, 

= e = e{Np + Ne) for E-B >0 and = e.f = ~e{Np + N^) for E-B <0. (7) 

The electric charge J° is given by the number densities of electrons and positrons as expected. On the other hand 
we need to explain why the form of the chirality is given as in Eq(I7]). Because positrons ( electrons ) are accelerated 
into the direction parallel ( anti-parallel ) to E, the momentum of positrons ( electrons ) is given hyp>0{p<0) 
foT E ■ B > OT p < { p > ) for E ■ B < 0; B — {0,0, B) with B > 0. { The particles are spontaneously produced 
with their momentum equal to because of the Pauli principle. ) Then, the helicity of both positrons and electrons 
is positive for £' • i? > or negative for £' • _B < 0, since the spins of the positrons ( electrons ) are pointed into the 
direction parallel ( anti-parallel ) to B. Therefore, the chirality is given as in Eq©. This is owing to the fact that 
the produced particles occupy the lowest Landau level with their magnetic moment parallel to B and their momenta 
are determined by the electric field. 
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III. CHIRAL ANOMALY AND PLASMA OSCILLATION 

Next, we explain the utility of the chiral anomaly and how the pair production under homogeneous electric and 
magnetic fields is described by using the anomaly. Suppose that the magnetic field is sufficiently strong for the 
charged particles to occupy only the lowest Landau level. Then, we may use the effectively two dimensional foumulas 
as described above. The anomaly equation is given in the following, 

dtJ'^^m-B (8) 

where we have taken an average over the transverse space xt = (xi, X2), of the chiral current so that / dxrdfj'^'^ = 0. 
We have also used the relation d^J^'^ = 0. The term in the right hand side of Eq® arises by taking account of 
quantum effects ( loop diagrams ) of electrons. That is, although the chiral current is conserved classically, the 
conservation is violated by the quantum effects. It is important to note that the violation term is given by the 
product of the electric and magnetic fields. 

When the electric field E is switched on at t = in vacuum, the pair production arises and the chiral charge 
is produced according to the anomaly equation ([8]). Suppose that the electric field is parallel to B. Then, positrons 
move to the direction parallel to E, while electrons move to the direction anti-parallel to E. We have Up cx 9{p) for 
positrons and rip cx 0{—p) for electrons. Thus, it follows from Eq(l5]) that = Np + Ne- Since Np — Ne in the pair 
production, the anomaly equation becomes 

dtN = sE (9) 

with N = Np and s = e'^B/Air'^. Similarly, when E is anti-parallel to B, then Up cx 0{—p) and hp cx 0{p). Thus, 
= —Np — A^e- Since E ■ B — —EB, the anomaly equation (|9]) is also obtained. Hence it holds in both cases with 
E ■ B ^ EB and E ■ B = -EB. 

Obviously, the number density of the electrons and positrons produced by the electric field is governed by the 
anomaly equation Q. For example, when E is switched on at i = 0, the particles are produced so that the number 
density is given by iV = sEt. This is the result with no back reaction of the charged particles on E. In order to take 
into account of the back reaction, we need to solve a Maxwell equation dtE ^ —J^ = —2eN as well as the anomaly 
equation the Maxwell equation holds in a spatially homogeneous system ( rotB = ). From these equations, we 
derive the equation 

d^Eit) + 2esE{t) - d^,E{t) + "^Eit) = 0, (10) 

whose solutions with the initial conditions E{t ^ 0) = Eq and N{t = 0) = are trivially obtained; E = Eq cos(-\/2es t) 
and N = Eq^J s/2es\n{\/2es t). The electric field shows a plasma oscillation in Fig.l with the frequency of the order 
of V eB ( ~ in the glasma ). Although we have solved only classical equations, the anomaly equation involves 
all of quantum effects associated with the pair productions. Thus, the chiral anomaly is very useful tool for the 
investigation of the Schwinger mechanism when both strong magnetic and electric field are present. 
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FIG. 1: electric field E{t) (dashed line) and number density of electrons N{t) (solid line) with arbitrary scale 
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We should make a comment on the scales involved in the discussion. The above result ( the plasma oscillation ) 
holds in the limit of eB ^ eE where only the states in the lowest Landau level are relevant. On the other hand when 
eB <C eE, the states in higher Landau levels become important since the energy |pF(i)| = \Jodt'eE\ of electrons 
can be much larger than V eB soon after their production. In the glasma we have eB ~ eE ^ ^ ™quark ( 
square of quark mass ). We argue that even in such a case the states in higher Landau levels are still not important 
according to the following reason. Most of the electrons are produced with their energies equal to 0. After their 
production, they are accelerated and their energies increase. The energies of the electrons produced at t = reach 

the energy |p_f| = /q dt'eE ~ eEQj^ ^ \J eB at the time of Qj^, the typical time scale of the oscillation. Then, 
they may make a transition to a state in a higher Landau level with the energy \feB. Such electrons are those which 
are produced just after the electric field is switched on. Most of electrons still have their energies less than V eB. So, 
the number of the electrons whose energies reach the critical energy V eB is still a fraction. Most of the electrons stay 
in the lowest Landau level. 

It has recently been shown [T6| that the vacuum decay by the pair production is mainly caused by massless particles 
in the lowest Landau level even if eE ~ eB. Although their analysis does not involve the back reaction of the particles 
produced, it supports our approximation that the relevant states for the decay of the electric field are those in the 
lowest Landau level. 



IV. PAIR PRODUCTION IN EXPANDING GLASMA 



Up to now we have discussed the pair production in the system where the electric and magnetic fields are infinitely 
extending in space. This does not correspond to the case of the real glasma, which are expanding after the heavy ion 
collisions. Thus we now proceed to show how E and N behave in a system where the fields and particles expand. 
Namely, we calculate their behaviors in t and 77 coordinates in which the background field E^ is uniform in r\. The 
behaviors of the fermions in the transverse directions specified by the coordinates x and y are identical to those in 
Cartesian coordinates. They are frozen in the states of the lowest Landau level. 

We first rewrite the Dirac equation, ^'^(idt — eA^) + "f^(idz — 6^3)^^ = by using the formulas, 

Oq — dardr + 9o?y9„ = coshrySr — ^- d„, dz = dzrdr + dz'qdr, = — sinYirjdr H ^ d„ (11) 

r r 

and 

sinh Tj cosh Tj 

^0 = Ar + dcti] A„ = coshrj Ar A„, Az — dzT A^ + dzTj A,, = — sinh 77 H A„. (12) 

T T 

Then, we obtain 

(j°{idr - eAr) + ^7^(i9^ - eA^)) </> (13) 

where we have set ij] = U(t>/y/T with U = cosh ^ + 7^7'^ sinh |. 

We solve the eauation p^ with Ar — Q hj assuming that A^i depends only on r. The assumption corresponds to 
the fact that the glasma is homogeneous in 77. The solutions are given by 



/ ■^^'^^P^*^'''^)'"^-^^ ^ith = exp(=F«fc,,logr =F«f^(T)) and = ' 

V / 



(14) 



where A^ are constants and J7 satisfies = eA^Jr. We have taken into account of the fact that (p or ^ represents 
only the states in the lowest Landau level. 

Using these solutions, we define the field operator corresponding to the field operator in eq([2|). 



'^^^ • exp 



2\/27r 



i2k^v)Ak{9ik,,)v+f+ + d{-k,-,)v-f-) + cM-^Kv)Bl{e{kr,)v+f\ + e{-kr,)v-fZ^) J (15) 
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where the first terms ( oc exp(jfc^ — i|fc^| logr) ) represent the solutions of the positive frequency and the second ones 
( oc exp(— ifc^ + «|/c,7| logr) ) do the solutions of the negative frequency. A), and satisfy the commutation relations, 

{Ak, AI,} — {Bk, bI,} = 2TrS{kn — k'^), the other commutators vanish. (16) 
Using the commutation relations, we find 

{UV, r)Al{v', r)} = Sa,bd{fj - 77'), (17) 

with a, = 1, or, 3. 

Now, we express the chiral anomaly in r and 77 coordinates. It is easy to see that 

doJ^'° + daJ^'^ = drJ^''' + - J^'^ + 9. JS'" = 2sr,E„ (18) 

T 

with J^'°- = {'ipj'^j^'ip) ( a = 0, 3 ) and Sr, = e^i?,,/47r^, where the last term drjJ^'"^ is assumed to vanish corresponding 
to the fact that the glasma is homogeneous in 77. The current J^'"^ can be represented in terms of J^'° and J^'^ in the 
following, 

J5.^ = doTJ^'° + darJ^'^ = cosh?7 J5'° - sinh?? J^'^ (19) 
Hereafter, we consider only the region of the central rapidity r/ ~ 0. Then, we have 

J5,r ^ j5,0 ^ (i^t^5^^_ (20) 
T 

In addition to the anomaly equation, we exploit a Maxwell equation, doF^'^ = J^, in order to take into account of 
the back reaction of the charged particles on Er,. The equation is represented in r and 77 coordinates as 

dArF^'^) = = drE, = -rJ^ ^ -J^ = -(^^Ij^), (21) 

where we have taken only the region of the central rapidity. 

Therefore, we have the following equations to find the temporal behaviors of the electric field E^ and the number 
density of positrons or electrons, 

drJ^'^ + -J^''' = 2sr,En and a^i; = -j3 = (22) 
T r 

with J^'"^ = (-^^7^(^), where we have assumed that the electric field is parallel to the magnetic field E-B = E,^B^ > 0. 
When the electric field is anti-parallel to the magnetic field, the anomaly equation becomes drJ^'^ + ^ J^'^ = — 2s^S^. 

In the above equations the expectation values of the currents J^'"^ and are taken by assuming appropriate states 
of positrons and electrons such that, 

J5.- = (i<^t^5^^ = / ^ (K - nk)e{kr,) - (nfc - hk)e{-kr,)) (23) 



') f dk \ 

-)= ^K^"'^ " nk)e{k^) - ink - fik)e{-krj)) , (24) 
where the state | ) is supposed to satisfy 

{AlAk') = 2nS{kr,-k'^)nk{t), {BlBk') = 2nS{kr, - k'^)hk{t) and {AIbI) = {BkAk') = 0. (25) 
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The quantities and ^ denote the momentum distributions of positrons and electrons, respectively. This can be 

easily understood by noting that the electric charge density is given by J° = ^"^"^^ = / ^^76(71^ — n/c). Obviously, 
J° = since electric charge is conserved and the initial state has no electrons and positrons. 

We should make a comment that when the electric field is parallel to the magnetic field, i.e. E ■ B ~ E,^B,^ > 0, 
positrons ( electrons ) in the pair production have momentum fc^ > ( fc^ < ) owing to the acceleration by the 

electric field. On the other hand, when E ■ B ~ —E^jB^f < 0, positrons have fc^ < and electrons do kjj > 0. That is, 

{uk ~ fik)e{k^) = nkd{k,,) ^ and {uk ~ fik)9{-k„) = -hk9{-k„) ^ for E ■ B > 0, (26) 
, {uk - nk)9{k^) = -fikOik^) ^ and {uk - fik)9{-k„) ^ nk9{~k^) ^ for E ■ B < 0. (27) 

Therefore, we find from these equations that the chirality J^'"^ is given by 

J^-'- = N + N for E-B>0 J^-^ = -{N + N) for E-B<0, (28) 

where N = 2iFt — ^^'^ ^ ~ 2irr — represent the number densities of positrons and electrons, respectively. 
Therefore, we obtain the following equations governing the pair productions, 

N 

drN + — = s^E,, and drE^ = -2eN (29) 

T 

where we have taken into account of the fact that the number of positrons is equal to that of electrons in the pair 
production, N = N. 

We can explicitly solve the equations ([29)) to obtain the temporal behaviors of the electric field i?,, and the number 
density of the charged particles. It is easy to see that Ejj satisfies 



d^E^ + -drE^ + 2es,,£;,, = 0, (30) 



and A'' is given by = —drE^i/^e. We obtain the solutions with the initial conditions, E.^{t — Q) = Eq and 
N{t) (r = 0) = 0, 



E = EoM^j2^r,T) and N= '' |Ji(72^t)| (31) 

where Jo,i denotes Bessel functions. We have shown the behaviors of these quantities in Fig. 2. We can see that they 
decrease oscillating with time t, while they simply oscillate without the decrease in the ordinary Cartesian coordinates 
( Fig.l ). The decrease comes from the expansion of the glasma. 




T 



FIG. 2: electric field E{t) (dashed line) and number density of electrons N[t) (solid line) with arbitrary scale 



The implication of the figure is in the following. When the electric field i?^ > is switched on, the pair production 
arises and the number density N increases, while the field becomes weak owing to the energy lose. The energy lose is 
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caused from the acceleration of the charged particles. There are two effects which make the number density increase 
or decrease. The number density increases by the pair production ( s^£'^ ), while the expansion ( —N/t ) makes the 
number density decrease: drN = SrjEjj — N/t. At the beginning of the pair production, N increases because the effect 
of the pair production is stronger than that of the expansion. But, at the time t = tq, N stops increasing because 
both effects balance with each other. After that, the effect of the expansion becomes stronger than that of the pair 
production. Thus, N begins to decrease. On the other hand, the electric field loses its energy with the acceleration of 
the particles and vanishes at the time r = ri > tq; E^i^t = ti) = 0. Then, the field changes its direction, i.e. < 0. 
The particles begin to be accelerated in the direction opposite to their velocity fc,, > 0. Thus, the pair annihilation 
begins and makes the number density decrease furthermore. Because the electric field accelerates the particles into 
the direction opposite to the particle velocity, the field gains the energy from the particles. Thus, the field becomes 
strong. Eventually the number TV of the particles vanishes ( N{t = r2) = ) at the time r = T2 > ri when the 
annihilation stops and the increase of the field strength also stops. We can see that at the time r = T2 a new pair 
production begins to take place under the effect of the electric field E,j{t = T2) < 0. This is the physical explanation 
of the behaviors depicted in the figure 2. 

We should also mention that the life time ti ~ 2.4/^2es^ ( the first zero point of the field E oc Jo(-y/2es^T) ) 
of the electric field is larger than the corrensponding one ti = 7r/2-\/2es of i? oc cos{\/2est) in the non-expanding 
glasma with s = Sr, ot B = B^i. Since the expansion makes the number density of the charged particles after their 
production lower than that in the non-expanding case, the rate of the energy loss of the field is slower than that in 
the non-expanding case. Hence, the life time ti is larger than ti. 

V. SUMMARY AND DISCUSSION 

To summarize, we have shown that due to the pair production of the charged particles, the electric field homogeneous 
in rj decreases with the oscillation in the r and rj coordinates. This should be contrasted with the case in the Cartesian 
coordinate where the field homogeneous in z simply oscillates without the decrease. The number density of the 
electrons shows the similar behavior as that of the electric field. Since r describes a proper time associated with the 
expanding fluid, it is natural that the energy of the electric field and the number density become lower with time r 
for 77 fixed. This lower number density of the charged particles leads to the longer life time of the electric field. We 
expect that these behaviors shown in QED also arise in the real glasma produced in high energy heavy ion collisions. 

We have only discuss the quark pair production under the color electric field produced in high energy heavy ion 
collisions. In the collision the most important products would be gluons, which can be also produced in the Schwinger 
mechanism[l^. In order to discuss the gluon production, we need to take into account of the effect of the color magnetic 
field. As is well known, the gluons are unstable under the magnetic field. It means that the gluons are coherently 
produced [20*1. Just as Higgs field 4>h located at the top of the potential 4>h = 0, the gluon fields exponentially grow 
to approach the bottom of the potential. The process can be regarded as the coherent production of the gluons. It 
is important to find the ratio of the number density of the gluon to that of the quarks produced in the Schwinger 
mechansm. In near future, we discuss the gluon production as well as the quark production by taking account of the 
back reaction of the coherent gluons and the quarks. 



We would like to express thanks to Dr. K. Fukushima for making me pay attention to the recent paper by Berges 
and Sexty. 
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